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We present a method for extracting tunnelling amplitudes from perturbation expansions which 
are always divergent and not Borel-summable. We show that they can be evaluated by an analytic 
continuation of variational perturbation theory. The power of the method is illustrated by calcu- 
lating the imaginary parts of the partition function of the anharmonic oscillator in zero spacetime 
dimensions and of the ground state energy of the anharmonic oscillator for all negative values of the 
coupling constant g and show that they are in excellent agreement with the exactly known values. 
As a highlight of the theory we recover from the divergent perturbation expansion of the tunnelling 
amplitude the action of the instanton and the effects of higher loop fluctuations around it. 



I. INTRODUCTION 

1. Tunnelling processes govern many important phys- 
ical phenomena. Their theoretical description requires 
the calculation of the contribution of critical bubbles to 
the partition function, including their fluctuation entropy 
0. The latter can be found at most to one-loop order. 
Higher loop effects are prohibitively complicated 0- It 
would be of great advantage tunnelling amplitudes could 
be derived from ordinary perturbation expansions around 
the free theory since these can be performed to many 
loops 0. The difficulty arising in such a program is that 
tunnelling amplitudes are described by the analytic con- 
tinuation of divergent Borel-summable power series ex- 
pansions in the coupling constant g to negative g where 
they become non-Borel-summable. None of the currently 
known resummation schemes 0, is able to deal with 
such expansions. Some time ago it was suggested that a 
resummation is possible by variational perturbation the- 
ory However, the imaginary parts calculated there 
gave accurate imaginary parts only in the sliding regime 
of larger negative g and did not invade into the proper 
tunnelling regime of small g dominated by critical bub- 
bles. A separate variational treatment of an instanton 
calculation was set up to cover this region ^6f . 

In this note we show that the non-Borel-summable se- 
ries can be evaluated with any desired precision by an 
appropriate continuation of variational perturbation the- 
ory @,|3j in to the complex g-plane. 

Variational perturbation theory has a long history 
S El EH • It is based on the introduction of a dummy 
variational parameter on which the full perturbation 
expansion does not depend, while the truncated expan- 
sion does. An optimal f2 is then selected by the principle 
of minimal sensitivity , requiring the quantity of inter- 
est to be stationary as function of the variational param- 
eter. The optimal £1 is usually taken from a zero of the 
derivative with respect to fl. If the first derivative has no 
zero, a zero of the second derivative is chosen. For Borel- 
summable series, these zeros are always real, in contrast 
to statements in the literature [13, Qj Il3 advocating 
the use of complex zeros. Complex zeros, however, pro- 



duce in general wrong results for Borel-summable series, 
as was recently demonstrated in Ref. 16]. 

The purpose of this note is to show that there does 
exist a wide range of applications of complex zeros in 
the previously untreatable field of non-Borel-summable 
series. These arise typically in tunnelling problems, and 
we shall see that variational perturbation theory provides 
us with an efficient method for evaluating these series and 
rendering their real and imaginary parts with any desir- 
able accuracy if only enough perturbation coefficients are 
available. The choice of the complex zeros is dictated by 
the requirement to achieve at each order the least oscillat- 
ing imaginary part when approaching the tip of the cut. 
We shall call this selection rule the principle of minimal 
sensitivity and oscillations. 

2. For an introduction to the method consider the ex- 
actly known partition function of an anharmonic oscilla- 
tor in zero spacetime dimensions, which is a simple inte- 
gral representation of a modified Bessel function K v (z): 

1 f°° 

Z(g) = -= / dx exp (-x 2 /2 - g x 4 /4) 

= e^i^g)- 1 ' 2 K 1/A (l/8g) , (1.1) 

For small g < 0, the function Z(g) and its inverse 
D(g) = Z~ x {g) have a divergent non-Borel-summable 
power series: 

oo 

D(g)=J2 a i9 l > (1-2) 

Z=0 

In the strong-coupling regime, there exists a convergent 
expansion 

oo 

D(g) = g<* J2 h 9~ U , (1-3) 
;=o 

with a = 1/4 and w = 1/2. In the context of critical 
phenomena, the exponent ui coincides with the the Weg- 
ner exponent 0] of approach to scaling [Isj . The Lth 
variational approximation depending on the variational 
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parameter Q is given by the truncated series 0, Q 

L 

C>(^)^ P E^)'^)' ( L4 ) 

where q = 2/w = 4, p — aq = 1. Introducing the param- 
eter cr = il 9 ~ 2 (f2 2 — l)/<7, the re-expansion coefficients 
are 



^-^/y-ay-K (i.5) 



Following the principle of minimal sensitivity, we have 
to find the zeros of the derivative of dnDyQ (g, Q), which 
happen to coincide with the zeros of a function of the 
variable a only, to be denoted by (a) : 

C (L V) = X> (^- lq) / 2 )(p-l q +2l-2L)(-a)^ 



1=0 



L-l 



(1.6) 



For a proof of this remarkable property see [n|, and the 
textbook (p. 291). 

As an example, we take the weak-coupling expansion 
to L = 16th order and calculate real and imaginary parts 
for the non-Borel region — 2 < g < —0.008 selecting the 
zero of Q le > (cr) according to the principle of minimal sen- 
sitivity and oscillations. The result is shown in Fig. ^ 
In order to point out how well the variational result ap- 
proximates the essential singularity of the imaginary part 
oc exp(— l/4g) for small g < 0, we have removed this fac- 
tor. The agreement is excellent down to very small —g. 

3. Let us now turn to the nontrivial problem of summing 
the instanton region of the anharmonic oscillator for g < 
0. The divergent weak-coupling perturbation expansion 
for the ground state energy of the anharmonic oscillator 
with the potential V(x) = x 2 /2 + gx A is, to order L; 



E, 



(L) 

0,wcak 



(<?) 



L 

E 

1=0 



ai g 



(1.7) 



where m = (1/2, 3/4, - 21/8, 333/16, 
30885/128, . . . ). The expansion is obviously not Borel- 
summable for g < 0, but will now be evaluated with 
our new technique, proceeding in the same way as for 
the above test function D{g) via Eqs. I|1.4fl through 
(|l.ti|) . The known strong-coupling growth parameters are 
a = 1/3 and u) — 2/3, so that p = 1 and q = 3 in 
Eq. l|1.6f) which will guarantee the correct scaling prop- 
erties for g — > oo. To order L = 64 we obtain from 
the optimal zero of C^ 64 - 1 (a) the logarithm of the scaled 
imaginary part 



1(g) := log V^572 E^ T (g) - l/3g, 



(1.8) 
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FIG. 1: Imaginary and real parts from variational perturba- 
tion theory of 16th order Dvir'(ff) as a function of log(— g) 
(dots) in comparison with exact curves (curves). In the imag- 
inary part we have removed the leading essential singularity 
by dividing out a factor exp(l/4<?) to see the amazing accu- 
racy with which this singularity is approximated. For very 
small g the onset of oscillations in the imaginary part can be 
seen which moves towards the origin for increasing order L of 
the variational approximation. 



shown in Fig. for -0.2 < g < -0.006. The point g = 
—0.006 is the closest approach to the singularity at g = 
for L — 64 before the onset of oscillations. 

Let us compare our curve with the expansion derived 
from instanton calculations 1211 : 



f(g) = h9 ~ W + b 3 g 3 - hg 4 



(1.9) 



with coefficients bi = 3.95833, b 2 = 19.344, b 3 = 
174.21, 64 = 2177. This expansion is divergent and non- 
Borel-summable for g < 0. Remarkably, we are able to 
extract this expansion from our data points. Since our 
result is given by a convergent expansion, the fitting pro- 
cedure will depend somewhat on the interval chosen over 
which we fit our curve by a power series. A compromise 
between a sufficiently long interval and the runaway of 
the divergent instanton expansion is obtained for a lower 



3 











• 










• 




• 




Kg) 



































-2 -3 -4 -5 

log(-ff) 



FIG. 2: Logarithm of the imaginary part of the ground 
state energy of the anharmonic oscillator with the essen- 
tial singularity factored out for better visualization, 1(g) = 

log \^/ —ng /2 S^varCsO] ™l/33> plotted against small negative 
values of the coupling constant —0.2 < g < —.006 where the 
series is Borel-nonsummable. The thin curve represents the 
divergent expansion around an instanton of Ref. |2ll| . The fat 
curve is the 22nd order approximation of the strong-coupling 
expansion, analytically continued to negative g in the sliding 
regime calculated in Chapter 17 of the textbook Q{. 



limit g > -.0229 ± .0003 and an upper limit g = -0.006. 



Fitting a polynomial to the data, we extract the following 
first three coefficients : 

bi =3.9586 ± .0003, 62 = 19.4 ±.12, 63 = 135 ± 18(1.10) 

4. The agreement of our curves with the exact results in 
Figs, n an d of our expansion coefficients in (|1.10|) with 
the exact ones in (ll.9fl demonstrates that our method is 
capable of probing deeply into the instanton region of the 
coupling constant. 

Let us end by remarking that another procedure of 
summing non-Borel series can be deduced from a devel- 
opment in the first of Refs. 01 ( see also Chapter 20 of the 
textbook 0]) One may derive a strong-coupling expan- 
sion of the type l|1.3fl from the divergent weak-coupling 
expansion, which can then be continued analytically to 
negative g by a simple rotation of the power g~ ul to 
e -vKul^_g^-wl _ xhis method is applicable only in the 
sliding regime. In Fig. [21 we have plotted the result- 
ing curve to order L = 9. The present work fills the 
gap between the sliding regime and the tunnelling regime 
by extending variational perturbation theory to all g ar- 
bitrarily close to zero, without the need for a separate 
treatment of the tunnelling regime. 

There exists, of course, a wealth of possible applica- 
tions of this simple theory, in particular to quantum field 
theory where variational perturbation theory has so far 
yielded the most accurate critical exponents from Borel- 
summable series 0, E, H3 • 
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